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Abstract. Energy transfer plays a vital role in many natural and technological 
processes. In this work, we study the effects of mechanical motion on the excitation 
transfer through a chain of interacting molecules with application to biological scenarios 
of transfer processes. Our investigation demonstrates that, for various types of 
mechanical oscillations, the transfer efficiency is significantly enhanced over that of 
comparable static configurations. This enhancement is a genuine quantum signature, 
and requires the collaborative interplay between the quantum-coherent evolution of 
the excitation and the mechanical motion of the molecules; it has no analogue in the 
classical incoherent energy transfer. This effect may not only occur naturally, but it 
could be exploited in artificially designed systems to optimize transport processes. As 
an application, we discuss a simple and hence robust control technique. 
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1. Introduction 

Transport phenomena can be safely enumerated amongst the most important physical 
processes occurring in both natural and artificial systems, and therefore the transport 
of charges or energy plays a central role in many scientific disciplines. Particularly 
important examples are the life-enabling transport processes in the molecular machinery 
of biological systems, which take place at scales ranging from a few atoms to large 
macro-molecular structures such as photosynthetic complexes. Pioneering studies of 
energy transfer, both theoretical and experimental, have already appeared in the early 
stages of the development of quantum mechanics [T], and this subject has been of 
interest ever since. Recent experimental investigations reveal the presence of long- 
lived quantum coherence during charge or energy transport processes in photosynthetic 
complexes [21 El IH E], charge transport through DNA [6], and in polymers [7j, even 
at physiological temperatures [H [9]. These findings raise the following question: To 
what extent may a fully quantum-coherent nature of the transport processes (rather 
than a classically incoherent, diffuse hopping), and perhaps even the presence of 
stronger quantum signatures such as quantum entanglement [TU] be responsible for the 
remarkable efficiencies that we witness in these systems [HI [121 lEl E] • 

At room temperature, motion of the underlying atomic and molecular structures, 
along which transport happens, is a ubiquitous phenomenon. Perhaps more important 
for living systems is the origin of molecular motion due to driving. The relevant time- 
scales range from a few femtoseconds for a stretching mode oscillation of two covalently 
bound atoms to several picoseconds for the collective modes of larger molecules, and 
beyond for the conformational changes of proteins [T^l [IE]- If the transport process 
and the underlying structure evolve on comparable time-scales, their dynamics do not 
separate, and we must expect an interplay between both that may affect the transport 
efficiency. 

In the present work, we treat an abstract transport model in which the motion 
of the molecule backbone structure is classical, and focus on the quantum-coherent 
transport that takes place on this moving structure. The quantum degrees of freedom 
that are carried by the molecular backbone will thus inherit time-dependent features 
from the classical motion. We show how such a time dependence can actively drive the 
coherent excitation transfer and thereby increase its efficiency beyond any comparable 
static scenario. The presence of this enhancement is a genuine quantum effect as it does 
not occur in a classical diffusive hopping transfer. With this we identify and quantify 
the motion- induced quantum enhancement of the quantum-coherent transport. 

It is known that dynamical structural changes in proteins affect transport processes 
in various ways. Environmental noise induced by thermal motion, and specifically 
dephasing, may lead to an enhancement of the efficiency of transport through a 
network of sites [T71 HH HH [201 ED |22]. There, the generic idea is that dephasing 
can overcome the trapping of energy due to a disorder of the local energies, or 
due to destructive interference along different paths through the network (similar 
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as destructive interference in a precisely aligned interferometer). As an example, 
such a dephasing process can result from coupling the local energies to the collective 
vibrations of the network. The transport process itself may, however, induce a 
rearrangement of the surrounding molecular structure that acts back on the transport 
process [211 ESI ESI EZl I2B] , and may induce a directionality, for example. In this 
case, protein motion and transport processes are coordinated as opposed to random 
changes due to thermal motion. Indeed, such triggered, long-lived coherent motion of 
protein nuclei has been found to exist [291 EO], even at room temperature. 

Here, we choose an approach consistent with the latter observation and treat 
the molecule motion as coordinated with the transfer process rather than as a noise 
source. Therefore, in our case the enhancement of transport efficiency has nothing 
to do with dephasing, but arises from driving, i.e. a constructive interplay between the 
structural changes of the transport network and the quantum-coherent hopping between 
the sites. Since the underling physical mechanisms of single-electron, hole, proton, 
or spin transport have an identical formal modeling, the results of this study can be 
straightforwardly applied to the corresponding parameter regimes of these systems as 
well. Furthermore, we demonstrate simple and therefore robust control techniques to 
induce or enhance these effects in linear systems. 

In the present work, we shall deal with a rather abstract model to investigate the 
effects of motion on transport efficiencies. Nevertheless, several biological applications 
are conceivable that largely resemble our chosen parameter regimes, and show features 
that agree with the assumptions underlying our model. 

In proteins, it is suggested that energy transport plays a crucial role for 
conformational changes [311 |32]. For the protein to change its shape, and perform 
its function, it is often necessary to overcome an energetic barrier. It is hypothesized 
that the necessary energy needs first to be available as vibrational energy in the protein, 
and that it is transported from a site where it is provided by ATP to the region within 
the protein where it is consumed by a conformational change. A type of secondary 
structure in proteins, in which the chain of amino-acids forms an a-helix, is suggested to 
provide such a transport pathway for vibrational energy [3ll|32]. Although the sites that 
carry the excitation in form of a localized molecular vibration (the amide-I oscillator is 
predominantly a C = stretching mode oscillation) are arranged in a helix, the structure 
is essentially linear. In this scenario, energy transport has been investigated under the 
name of the Davydov-Scott-model [211 ESI ISl] , to a greater degree of realism than we 
aim to achieve here, and with an emphasis on localization phenomena of excitations on 
the a-helix due to motion [SSI ESI Elj ■ 

For photosynthesis in green sulfur bacteria, light is harvested by a large antenna 
complex, and the excitation is subsequently transferred via the well-studied Fenna- 
Matthews-Olson (FMO) protein, that thus functions as a wire, to the reaction center 
where the excitation energy is converted to chemical energy, for a review see [38] 
and references therein. Within a single unit of the FMO trimer, the excitation has 
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been shown to propagate quantum-coherently via seven bacteriochlorophyll pigments, 
entering and exiting at specific pigments [21 |3l O [8]. In the reaction center, an 
electron is then quantum-coherently transferred via several sites. This electron transfer 
is accompanied by a structural changes in form of low frequency vibrations of the 
surrounding protein environment on the same time scale [221 EOl S] • 

The paper is organized as follows. In section |2| we introduce the formal model, and 
start our investigation with the simplest possible case of energy transfer, namely the 
dimer model that consists of only two molecules. This serves as a reference system to 
build intuition of how mechanical oscillations can enhance the energy transfer. By 
comparing quantum excitation transfer with the classical Forster theory of energy 
transfer, we explain how mechanical motion can constructively cooperate with the 
quantum dynamics of the excitation. In section |3| we generalize this model to moving 
systems of many molecules, where we investigate motion modeled by their normal modes 
and the resulting effects on excitation transfer. Section |4] extends the study of motion of 
a chain towards a scheme in which the interaction strength between sites is modulated 
by a pulse. 

2. Transport Model 

Let us first set the formalism to describe excitation transport. Here, we limit ourselves 
to a linear chain of sites, which represents an array of molecules, e.g. an a-helix, 
among which an excitation can be exchanged due to dipole-dipole couplings between 
the molecules. When the probability of an excitation being present in the system is 
low, we can restrict our analysis to the single-excitation subspace. Furthermore, we 
assume for simplicity that the molecules are sufficiently distant from each other such that 
their interaction can be reduced to the dominant interaction between nearest neighbors. 
Generally, however, this is not the necessarily the case, as in the FMO complex, for 
example. With these assumptions, the Hamiltonian of the system can be written as 
the tight-binding Hamiltonian of an interacting A^-body system in the single-excitation 
manifold {h= 1), 

N N-1 

H = ^en\n){n\ + ^ Jn{\n){n + 1| + |n+ (1) 

n=l n=l 

where |n) denotes the state with the excitation at the n-th site having energy e„, and 
Jn is the coupling strength between the n-th and {n + l)-th molecule, as depicted in 
Figure |l| (left). 

We assume the initial state of the system at time t = to be a single excitation 
localized at the first site (left end of the chain in figure [T]), i.e. p(0) = This is 

in accordance with the idea that the excitation enters the transport network at specific 
sites as verified in the FMO complex |39j, or provided locally by ATP in the a-helix 
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Figure 1. Left: Schematic transport model of next-neighbor coupled chain with local 
dissipation and the sink attached to the last site. Right: Schematic energy diagram 
of two electronic states before and after excitation showing how excitation can induce 
motion of the nuclei that is coordinated with the launch of the excitation propagation 
in the complex [25] . 



scenario. For spectroscopic investigations, fs-pulses of the right wave length are able to 
excite specific sites with in a complex of chromophores as done in the FMO kind |2]. 

The motion of molecules will impose a time dependence onto the Hamiltonian ([T|. 
We treat the spatial positions of the chain constituents as classical quantities that follow 
well-defined trajectories [IT]. This semi-quantal approximation holds if the involved 
molecules are too large to show their quantum behavior, i.e. if the uncertainty with which 
the position of a molecule can be determined due to the classical thermal fluctuations 
is larger than the quantum width of the associated wave function in the position 
representation (see also |1Q]). For a molecule of mass M that is attached to a spring 
with spring constant K, and that is in contact with an environment at temperature T, 
the uncertainty of its position due to thermal fluctuations is Ax^^ = 2kBT/K, while 
the width of the (ground state) wave function is Axg„ = 2h/ \/ MK, with being 
the Boltzmann constant, and h the Planck constant. Therefore, Axth > Axqu if the 
temperature is high enough, i.e. if ksT > h^K/M, and we can focus on classical 
motion of the excitation carriers. For example, the typical critical temperature for an 
a-helix structure is 60 K (as obtained for the mass and spring constant in the effective 
one-dimensional model: i^T ~ 39 — 58.5 N/m, M ~ 5.7 x 10~^^ kg [33]) which is certainly 
much below room temperature. 

Motion of the individual molecules changes their relative distance in time and 
thereby modulates the distance-dependent dipolar coupling. At the same time, 
deformation of the molecules leads to a change in dipole moments and thus also induces 
a time-dependence of the coupling. The latter will be dominant whenever the molecules 
are tightly embedded in a protein scaffold such as in the FMO complex. For our 
purposes, and in order to demonstrate the effect, we only focus on a modulation of 
Jn due to the change in distance, since we expect a change of dipoles in strength or 
direction to yield similar results. For concreteness, we choose a simple form of motion 
where neighboring molecules change their distance periodically, which holds for small 
amplitudes in the harmonic regime. We will later see that for typical times that the 
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excitation spends in the complex only few molecule oscillation periods at the relevant 
frequencies suffice such that we essentially disregard damping and dephasing of the 
molecule motion. For times short compared to the damping of the motion, the distance 
between site n and n + 1 is then given by 

dn{t) =do- [un{t) - M„+i(t)] = do[l- 2a„ sin(wt + , (2) 

where Un{t) is the displacement of the n-th molecule, do is the equilibrium distance 
between two neighboring sites, and a„ is the individual sites' relative amplitude of 
oscillation when they move with opposite phase around their equilibrium position. This 
affects the dipolar coupling strengths J„ according to 

7 (t) = = ^5 (s) 

[dnitW [l-2a„sin(u;t + 0„)]3' 

where Jq contains the dipole moments (here assumed to be constant) and physical 
constants, and we define Jo = Jo/d^, which has the unit of an energy. The extrema 
of the coupling are Jn,mm when the molecules are farthest apart from each other, that 
is, when dn = do{l + 2an), and Jn,max when their distance reaches the minimum value 
of dn = do(l — 2a„). The deterministic motion of the molecules imposed by ^ is 
for a given amplitude, frequency, and phase synchronized with the propagation of the 
excitation. This is in contrast to a stochastic motion of the molecules as caused by 
thermal fluctuations, for example. It is, however, conceivable that upon the arrival on 
the excitation in the complex the electron configuration changes and nuclei begin to move 
in this new potential due to the Franck- Condon principle as shown in figure [l] (right) [25] . 
The wave packet that describes the nuclear motion has been observed to exhibit 
surprisingly long coherence times in reaction center proteins [29], which is consistent 
with our ansatz to treat positions of the sites classically, and in concert with the initial 
excitation. 

Other infiuences that motion might have on the dynamics of the transport include 
the detuning of the excitation energies of the individual sites. Here, we will mostly 
ignore the such arising disorder and localization effects. In a realistic biological context 
such as excitation transport in light harvesting complexes, these effect are usually not 
neghgible [H]. 



In addition to the Hamiltonian evolution, the excitation dynamics takes place in an 
open quantum system and therefore also suffers from dissipation. We describe the loss 
of the energy excitation due to dissipation in Markov approximation by the following 
Lindblad super-operator: 

N 

Ldissp = 7n.[2a-p(T,+ - {a+(x-,p}]. (4) 

n=l 

where (a^) is the creation (annihilation) operator of the excitation at site n, and 
{A, B} = AB + BA. Since realistically each molecule feels a different environment. 
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the local dissipation rates 7„ are different in general. For simplicity and due to so far 
unavailable experimental data, however, we will later assume them to be equal. 

In order to measure how much of the excitation energy is transferred along the chain 
(and not lost due to dissipation), we introduce an additional site, the sink, representing 
the final (iV + l)-th trapping site that resembles a reaction center, for example. The sink 
is populated via irreversible decay of excitation from the last site, N . This approach for 
quantifying the energy transfer efficiency follows Ref. [T3|, and is formally implemented 
by adding the Lindblad operator 

to the master equation, where 75 denotes absorption rate of the sink. 

In order to calculate the efficiency of the energy transfer given by the asymptotic 
population of the sink, it is necessary to integrate the following master equation {h= 1): 
dp 

— = Lp = i[p, H{t)] + LdissP + -^^sinkP- (6) 

Since we are interested in the asymptotic population of the sink, which is our figure of 
merit for the transport efficiency, we search for long-time solutions of the above equation. 
Asymptotically in time, the excitation will be either lost to the environment due to L^\ss 
or, ideally, trapped in the sink due to Lsink- Note that it therefore does not suffice to 
merely search for the stationary state, which is characterized by the ground state of the 
sites: it thus contains no information about how likely the excitation is trapped in 
the sink instead of being lost to the environment. 

In numerical simulations, we will usually consider uniform local energies e„ = e 
unless otherwise noted. Additional effects that arise from static disorder or detuning of 
the local energies due to the coupling to the nuclear motion are outside the focus of the 
present work, and are thereby mostly neglected. Henceforth, all energies, time-scales, 
and rates will be expressed in units of Jq, and thereby we effectively set Jq = 1. 



Dimer Model 

In the remainder of this section, we study the case of a short chain composed of only 
N = 2 interacting sites, the dimer, plus the sink in a local dissipative environment for 
El = 62 = e. Although simple, this toy model exhibits all the essential features of larger 
systems, and gives us the possibility of guiding our intuition towards the mechanisms 
that underlie the transport processes. 

For a time-independent coupling J, the excitation coherently oscillates between 
both sites at a constant Rabi frequency until it is either lost into the environment or 
absorbed into the sink, when it is at site 2. The asymptotic sink population is given, in 
this static situation, by [19j 

^'^''^ ' (27 + 75)[7(7 + 75) + J2] ^ > 

where 71 = 72 = 7 is the rate of dissipation into the environment at each site, and the 
superscript "st" indicates the static, i.e. non-oscillatory, case with constant J. 




Figure 2. (a) Coupling strengtli. (b) Population dynamics for coherent transfer in a 
dimer molecule oscillating at the optimal frequency. Curves represent the population 
of site 2 (solid) , the sum of the populations in site 1 and 2 (dashed) , and the population 
of the sink (dotted). The time axis is rescaled using the oscillation frequency cj — 4.54, 
which corresponds to the maximum in figure pi a). 



In the case of a time- dependent coupling J{t), the master equation can be 
straightforwardly integrated if the combined dissipation rates at both sites are equal, 
i.e. for 7i = 72 + 75 = F, and ei = as above. With this choice, the populations of the 
respective sites read: 

Pi(t) = cos^ [J^ J{t')dt'^ e-^r*, and (8) 

P2{t) = sin2 J{t')dt'^ e-2r*. (9) 

The sink population is finally obtained by integration: Psink('^) = 275 P2{t)dt. 

Since the phase of the Rabi oscillation is obtained by integration over the coupling 
strength in ([s]) and we expect a faster propagation of the excitation between sites 
whenever the coupling strength is stronger, i.e. when the two sites are closer to one 
another. Similarly, we witness a slower excitation propagation (or even an almost static 
situation) during times of weak coupling, when the sites are distant. For a relative 
amplitude of a = ai = 1/4, figure |2](a) shows the time-dependence of the interaction 
strength, which exhibits (for the chosen a) short times of strong interaction, separated 
by longer durations of weak coupling. This asymmetry in the coupling strength is due 
to the disproportional dependence of J{t) on the relative distance of the sites, cf. 

For the right ratio of Rabi frequency and mechanical oscillation frequency u, the 
narrow peaks of strong interaction, when the sites are close, can be used to quickly 
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transfer the excitation between the sites, whereas the longer "valleys" of low interaction 
strengths and slow Rabi oscillation are used to effectively lock the excitation on one 
site, cf. figure [2|b). Using this quick-transfer-and- locking strategy, we can estimate 
parameters that allow us to maximize the duration during which the excitation is located 
at the second site, and thereby most efficiently exposed to the decay into the sink. 
Starting with the initial spatial position given by the phase (p = (pi = n /2 such that 
the sites are closest, and the interaction is strongest, we need about the first l/4th 
of the mechanical oscillation cycle, when the interaction is still sufficiently strong, to 
transfer the excitation from site 1 to site 2 (in figure [2]^b) P2 approaches the total 
population of both sides). During the next half cycle of the mechanical oscillation, 
when the interaction strength is low, the site population does not change significantly. 
This effectively amounts to a locking of the excitation at site 2. The following half-cycle, 
during which the interaction is strongest again, then needs to last long enough for the 
excitation to move again to site 1 and back. 

Comparing with ([s]) and ([9]), the above translates into the requirement that 
J^^^ J{t')dt' ~ 7r/2 for the first transfer to site 2, where T = 2'7i/u is the period of 
the mechanical oscillation. For a strong interaction modulation as in figure [2]^a), we 
neglect the coupling strength when the two sites are distant, and the integral can 
be approximated by Javg^r/w, where Javg = Jo(l + 2a^)/(l — 4a^)^/^ denotes the time- 
averaged coupling strength. For the above situation, we thereby obtain an optimal 
frequency of w ~ 2Javg. Higher-order maxima are obtained by adding full Rabi cycles 
to a transfer step, i.e. 

/ J(t')dt' ~ - + mvr and a; ~ m = 0,l,2,.... (10) 

Jo 2 2m + 1 

In contrast, for a locking of the excitation at site 1, i.e. for a minimum exposure to the 

sink decay, an analogous argument yields the minima of the transfer efficiency: 

J(t')dt'~m7r and w ~ m = l,2,.... (11) 

m 

Figure [3|^a) summarizes the transfer efficiencies for a range of mechanical oscillation 
frequencies. The estimated frequencies for extremal transfer agree well with the obtained 
numerical integration, even for the situation in the plot, in which the decay rates do not 
fulfill the assumption underlying the analytical solution in ([s]) and 

In case the combined decay rates are not equal at all the sites, the symmetry 
of the problem is broken, and the eigenstates of the Hamiltonian, the symmetric and 
anti-symmetric superposition of the excitation being at sites 1 and 2, are coupled to 
one-another by the decay process. Moreover, if the initial phase is not tt/2, maxima 
of Psink(i^) are no longer characterized by a perfect locking of the excitation at site 2 
for all times, but only during the first few relevant oscillations until the excitation has 
essentially left the system. 



In order to quantitatively demonstrate how the presence of mechanical oscillations 
genuinely enhances the energy transfer, we compare the asymptotic population of the 
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Figure 3. Asymptotic sink; population dependent on the mechanical oscillation 
frequency uj. (a) Initial phase of the mechanical oscillation is = tt/2, i.e. the 
first site of the chain is excited when the sites are closest. Estimated frequencies 



according to (10) and (11) for extremal sink population are indicated by vertical 
lines. Horizontal lines indicate the asymptotic sink population for static couplings. 
For asymptotically large w, the curve approaches the sink population Pginki'^'^vg) ■ 
(b) Amplitude dependence of the maximal enhancement with an initial phase (j) — ti /2 
(solid) and the optimal initial phase (dashed) shows an algebraic decrease for small 
amplitudes in a double logarithmic plot. The cross marks the amplitude used for the 
other plots. 



sink Psink(i^) for a chain oscillating with frequency u with the obtained sink population 
for a fixed chain in the closest and hence strongest coupled site configuration, -Pgfiik(t/max)- 
For this purpose we introduce the motion-induced quantum enhancement 

A(a;) = Psink(t^) - i^lk(^max). (12) 

We choose a set of parameters for all the following simulations that allow us to 
conveniently identify the effects. The chosen time scales are approximately of the same 
order: Jq = 1, 7„ = 1/10, and set 75 = 1/2. Figure [Sj^a) shows that for the moving 
molecules the sink population Psink(w) surpasses the maximal static case -Ps^nk( Jmax) by 
several percent for the chosen parameters, and exhibits an enhancement A > 0. 

We emphasize that the motion-induced enhancement beyond the efficiency of the 
static case at maximal coupling strength is a true quantum effect. That is, the motion- 
induced efficiency enhancement in the quantum mechanical case is stronger than the 
motion-induced enhancement in the classical case. The coherent excitation transfer 
(Rabi oscillation) can transfer all of the excitation to the second site. A modulation of 
the coupling strength on a comparable time-scale may then prolong the effective time 
that the excitation is localized close to the sink as compared to the static chain. In the 
corresponding classical case, however, the excitation transfer takes place according to 
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Forster theory |12], and is governed by a rate equation, e.g. 




= ^M„„(t)P„(t), 



(13) 



m=l 



where the population of the ra-th site P„(t) depends on the transfer rates Mnm for a 
stochastic hopping between sites n and m. For detailed balance Mnm = ^mn, this leads 
to a diffusive propagation of the excitation through the network of sites, and in absence 
of dissipation to an equilibration of the site populations. The maximal population at the 
last site that is connected to the sink is therefore at best P/v(t) = A modulation of 
hopping rates Mnmif) does only affect the speed with which this equilibration is reached 
but not the maximal population of the last site. The mechanical oscillation of the sites 
which modulates the hopping rates can therefore not achieve a hopping rate above the 
closest site configuration, where the maximal population of the last site is reached fastest, 
and maintained highest in presence of dissipation and decay to the sink. The nature 
of diffusive transport does therefore not allow for dynamic excitation locking strategies 
as pictured earlier. The corresponding quantity Aci(a;) = -Psink,ci(i^) ~ -fslnk ci('^max) 
then measures the efficiency of classical diffusive excitation transport through a moving 
molecule structure in comparison to classical transport through a static, maximally 
coupled structure. We can therefore conclude that a value of A(ci)(a;) > is classically 
forbidden. A positive motion-induced quantum enhancement A(a;) > thus proves and 
quantifies a "quantum advantage" due to the quantum nature of the transport process. 

The observed motion-induced efficiency enhancement is present for all values of 
amplitudes a and magnitudes of the modulation of the coupling J. However, the effect 
is less pronounced for smaller values than the ones used here. The optimum oscillating 
frequency u is then no longer characterized by a complete locking of the excitation 
at the second site, but rather by a modulation of the Rabi oscillation such that the 
excitation lingers slightly longer at the last site. Our particular choice of a = 1/4 
amplifies the visibility of the effect, but it does not cause its presence. Figure |3](b) 
collects the maximal enhancement A(a;opt) at the optimal frequency for several values 
of the amplitude a. Even for reasonably small amplitudes of the displacement, the 
effect is of the order of a percent. The doubly logarithmic plot suggests that the effect 
diminishes only algebraically for amplitudes ranging from large displacements that are 
expected for a driven motion, down to small position fluctuations as caused by thermal 
effects. 

In a-helices the equilibrium distance between the amide-I oscillators is do ~ 4.5 A 
whereas their relative displacement reaches values of doa ~ 0.1 A [331 |3T], which amounts 
to an a of one order of magnitude smaller. Nevertheless, the coupling J exhibits a 
variation of ~ 7%. For smaller amplitudes of the mechanical oscillation than the ones 
chosen to obtain figure [3| the extrema persist but are less pronounced. It is therefore 
conceivable that, at the right frequency of motion, even in biological or chemical systems 
where oscillation amplitudes are small, the quantum-coherent nature of this transport 
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Figure 4. Asymptotic sink; population Psink('^) for different initial phases of the 
mechanical oscillation. Thin grey curves correspond to 50 initial phases uniformly 
spaced in g [0, 27r]. Thick black curves indicate the enveloping curves (outer curves), 
and the phase averaged asymptotic sink population (inner curve) of 200 uniformly 
sampled phases. Horizontal lines again show the sink population for static cases 
with Jmax, ^avg, and Jo (from top to bottom). For asymptotically large lu, all curves 
approach the sink population of the static case with Javg- 




process on a moving structure yields efficiency improvements beyond what is classically 
achievable. Please also note that the total efficiency gain with respect to the equilibrium 
situation is much larger. 

So far, the discussion was limited to the motion of the molecules starting with the 
initial phase = 7r/2, i.e. with the molecule in the closest and hence maximally coupled 
configuration. This is justified - as explained earlier and shown in figure [l] (right) - 
when the motion is triggered by the excitation of the complex that changes the nuclear 
equilibrium position to greater distances than in the ground state. For motion of mostly 
thermal origin, such a coordination is generally not the case, and the initial phase 
would be random. In particular, one cannot expect an enhancement if the molecules 
are initially in a distant configuration and move away from each other. We therefore 
expect that, in the absence of synchronization between excitation transfer and motion, 
the molecules will not be able to efficiently drive the transfer. 

In order to test whether or not an enhancement of the excitation transfer above 
any comparable static case persists for an arbitrary initial phase of the mechanical 
oscillation, we statistically sample uniformly over initial phases. In figure |4], each thin 
gray curve represents the sink population and enhancement for a given initial phase of 
the oscillation, of which the figure |3] represents an instance. The outer black curves 
are the envelopes of the data, and represent the maximal and minimal sink population 
that can be obtained for the chosen parameters at a given oscillation frequency u. The 
inner black curve shows the average sink population over all initial phases. We can 
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observe that the average sink population for a random initial phase lies below the static 
case with coupling Jo for low frequencies, and approaches the efficiency of the static 
case with the time-averaged coupling strength Javg from below for large frequencies. 
An enhancement of the transfer efficiency over the static case with Jq can on average 
still exist simply due to the fact that the average coupling strength is stronger than Jq. 
Furthermore, it becomes evident that an enhancement A > can only be reached in a 
window of frequencies. For very slow motion, i.e. for small cj, the excitation propagates 
on a quasi-static structure, and a certain threshold of the underlying motion has to be 
passed before A > appears. For very fast motion, w ^ Jq the transport efficiencies 
converge towards the static configuration with the time-averaged coupling strength and 
thus A < (not shown in figure |4]). We can thus conclude that for a given oscillation 
frequency uj ^ 0{Jq), a large enhancement over the static situation can be obtained, 
but strongly depends on the initial phase. Therefore, it is crucial that the transport 
process can be initialized with a specific range of phases of the site oscillation rather 
than a uniformly distributed random phase. At the same time, however, if one has an 
experimental way of setting or infiuencing the initial phase, the strong phase-dependence 
might be exploited to switch or to observe this effect experimentally. 

In the following sections, we consider longer chains and perform a deeper analysis 
of the behavior of A(u;). We anticipate that the main features identified in the dimer 
model are qualitatively present also in larger systems, and that the quick-transfer-and- 
locking strategy can be employed to shuttle the excitation along a longer structure. In 
particular, we expect that A(a;) exhibits an analogous oscillatory behavior that extends 
to positive values for a specific range of frequencies. The identification of such optimal 
frequencies could in principle lead to a better control and to an enhancement of the 
efficiency of energy transfer in artificially designed systems via the modulation of the 
oscillation frequency. 

3. Multi-molecule oscillations 

In this section, we consider longer linear chains as a first step towards a more suitable 
representation of nature-inspired systems, such as the amide-I oscillators along an a- 
helix or the chromophores in a light harvesting complex. The analysis performed in 
the dimer case is now extended to linear chains of molecules. The underlying 
intuition is that the molecules carrying the excitation are embedded in a larger complex, 
e.g. a protein, but again are assumed to follow a concerted motion. The particular 
environmental and boundary conditions change from instance to instance. For this 
reason and for simplicity, we model the linear chain as a succession of molecules having 
the same mass and attached via a spring to their respective neighbors, as e.g. inspired 
by the amide units in an a-helix. For the first and last site we consider two possible 
situations: both the ends are (i) completely free or (ii) both are in contact with two 
"walls" via two additional springs. In the first case, the chain would represent an a- 




Figure 5. Left: Sloshing mode of oscillation in which half of the system of two-level 
sites contracts and half expands alternately. Right: Breathing mode in which the entire 
system contracts and expands periodically. The individual sites are coherently coupled, 
and the sink is connected via irreversible decay to the last site. 



helix structure without a narrow confinement in the longitudinal direction, while in the 
latter case the form of the surrounding protein determines an effective space in which 
the helix is confined. For both situations, we study collective normal mode oscillations 
of low frequency as a relevant example of naturally occurring oscillations. 

The transport efficiencies along the chain that are obtained under these dynamical 
conditions are compared to the ones in the static case, in which all the molecules 
are at the equilibrium distance. We thus obtain the gain that is simply due to the 
presence of motion. Alternatively, we compare to a configuration where the sites are 
maximally coupled, which yields A. In the case of uniformly coupled, static chains 
of more than three sites, the excitation is no longer completely transferred to the last 
site with probability one, from where it may be trapped into the sink [13]. Since, in 
addition, with increasing chain length the excitation is also more likely to decay into the 
environment, we expect lower absolute efficiencies than in the dimer toy model. On the 
other hand, there is the possibility that the motion-induced enhancement might add up 
over several transfer step along the chain. 

Before comparing the asymptotic amounts of population in the sink, we briefly 
introduce the model that we use to derive the proper motion of the linear chains. For 
simplicity, we assume that the forces between two neighboring sites are proportional 
to their relative displacements from the equilibrium positions: this approximation, 
known as the harmonic approximation, holds as long as the displacements are relatively 
small. To help in visualizing the system, one may think of the molecules forming the 
unidimensional lattice as connected by elastic springs (see figure [s]). In the chain of 
coupled oscillators, the internal force exerted on the n-th site is: 

Fn = M-^ = K{Un^i - Un) + K{Un+l - (14) 

where K is the spring constant (intermolecular force), M is the particle mass, and Un 
the displacement of the n-th site from its equilibrium position. For the first and last 
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site, we alternatively have only one spring or fix them with additional springs to the 
confining boundaries. At this point we neglect the infiuence of random thermal forces on 
the particles. An additional noise term in the position and thus in the coupling strength, 
as long as it is sufficiently small, will only slightly affect the effective Rabi frequency as 
compared to the case without noise, since the Rabi frequency it is effectively obtained by 
integrating over the coupling strength. The particular shape of J„(t), and fluctuations of 
it will therefore not greatly influence the transport efficiency as long as a quick-transfer- 
and-locking strategy can be achieved. 

A chain of coupled oscillators has a generic solution of the above equations of 
motion that can be expressed as a linear combination of normal modes. The normal 
modes are independent, collective modes of oscillation in which all the sites move with 
the same periodicity and do not cause the excitation of other oscillatory motions. Each of 
these modes is characterized by a frequency and a relative phase between the individual 
sites. For the confined chain, the normal mode frequencies and the displacement of the 
n-th site are (see e.g. [S]) 
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respectively, where cjq = \/K/M ^ Ag = sin[g7r/(A^ + 1)] is a normalization constant, and 
do is the equilibrium distance between neighboring sites. For the chain with open ends 
the expressions are 
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with Bq = cos[g7r/(2A^)]. Figure [s] illustrates the normal mode with the lowest frequency 
for the confined chain, the sloshing mode, and for the open chain, the stretching or 
breathing mode. 



For a more detailed study, we first analyze the effect of oscillations for individual 
normal modes of the confined chain. This approach is consistent with the fact that 
slow coherent motion of the underlying structure is indeed observed in proteins where 
transport phenomena occur [2^ [50] . We mimic this slow coherent motion by the low 
frequency collective modes of the chain. The sloshing mode, the lowest lying of the 
normal modes, represents a contraction of one half of the system, and a simultaneous 
expansion of the other half. For higher modes, the system is divided into more parts, 
each of which can be considered to be in a local sloshing mode. Figure |6] shows the 
increase of the asymptotic sink population over the static case for selected normal modes 
of a confined chain of coupled oscillators, when varying uq, i.e. spring constant and/or 
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Figure 6. Motion induced energy transfer gain over the uniformly coupled chain for 
different normal modes of TV = 13 confined sites. Oscillation amplitude constant is 
ttq = 1/24 and 4>q — tt. Only the modes q = 1,2, and 13 are shown. Jo) ~ 0.17, 

and the flat part of g = 13 for low wq amounts to a vanishing absolute sink population. 



mass of the molecules. The structural features that we found for the simple dimer case 
persist, and exhibit frequencies with increased and suppressed transfer efficiency. 

For low frequencies, a suppression in the transport efficiency dominates the 
behavior. In this frequency regime, the mechanical oscillation is almost static compared 
to the timescale of the coherent excitation transfer. The excitation transfer thus happens 
on a chain with disordered couplings, which lead to localization of the excitation at the 
beginning of the chain, and hence to a decreased transfer efficiency. 

For higher frequencies uo, we obtain efficiency maxima that surpass the transport 
efficiency of the static, ordered chain. As in the dimer model, these frequencies are in 
resonance with the wave-like propagation of the excitation through the chain. 

For asymptotically large uo all curves converge to the efficiency of the static case 
with the time-averaged coupling strength between nearest neighbors (not shown in 
figure [6]). Since we fix a, = 1/24 for all modes, and the relative distance between two 
sites is also determined by a position- and mode-dependent prefactor in (16) and (18), 
the time-averaged coupling strengths differ between the normal modes, and between 
different sites of a single mode. 



In order to elucidate the dependence of the transport efficiency on the parameters that 
define the system, we investigate different amplitudes and chain lengths for a single 
normal mode. Since the sloshing mode has no configuration where all the sites are 
closest and hence maximally coupled (half of the chain expands whereas the other 
half is compressed), we choose the breathing mode of the unconfined chain, where 
there exists such a configuration. For the breathing mode, figure [7] collects the data 
of the maximum enhancement of the transfer efficiency over the maximally coupled 
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Figure 7. Left: Maximum energy transfer enhancement over the maximally coupled 
static chain Aopt (top) and the corresponding optimal oscillation frequency cjopt 
(bottom) for the breathing mode for different chain lengths N and oscillation amplitude 
constants a = ai. Right: Frequency dependence (w = wi) of the energy transfer 
enhancement as compared to the chain at rest, i.e. for Jq, with amplitude constants 
Oi — N/12 (top), fli — 7V/48 (bottom), and (t>i = 0. The dotted lines include a 
detuning of the site energies due to the motion with coupling strength x = 10- 

static case, i.e. Aopt = max^-^ Psink(i^i) — -Psink(^n,max), and the corresponding optimal 
oscillation frequency Wopt- The amplitude constant a = ai is chosen such that the 
maximal extension of the chain increases linearly with A^, and the relative change in 
length of the entire chain during the mechanical oscillation, AL/L, remains constant 
when changing N . In particular, the case of a = provides an instance where the 

relative distance between neighboring molecules is of the order of what is also found in 
vibrating a-helices at room temperature . 

As for the dimer molecule, a positive Aopt cannot be achieved in the analogous 
classical description, and likewise here, for the quantum-coherent transport, positive 
values exist over a large range of chain lengths. We observe that, in general, Aopt 
decreases with increasing N . This is essentially due to the larger distance that the 
excitation needs to travel, and hence a proportionally increasing probability that the 
excitation is dissipated. We attribute the increase for short chain lengths to the 
accumulation of the efficiency gain, which is observed for two sites, over several sites. 
The decrease of Wopt with increasing chain length qualitatively follows the inverse of the 
time at which the excitation is (partially) transferred to the last site in a homogeneously 
coupled chain, i.e. with increasing chain length, the time also increases at which the first 
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Figure 8. Left: Decay of the motion-induced transport enhancement A for the 
breathing mode as a function of the dephasing rate for different chain lengths. 
Amplitude constant is a = ai = iV/48, (pi = 0, at a fixed frequency wq that is optimal 
at 7 = 0. Right: Critical dephasing rate 7c dependent on the chain length N . 



maximum of the last site's population appears [IS]. We can apply a similarly intuitive 
argument to motivate the dependency of cjopt on the oscillation amplitude factor a. 
When the amplitude factor a of the oscillation is large, the inter-site coupling strength 
is subject to a large variation and therefore, when the molecules are closest to each other, 
the excitation propagates very fast such that a quick modulation of coupling strength 
is necessary to prevent the immediate reversion of the population. 

We expect that this non-classical enhancement induced by oscillations can be 
observed not only for the chains with uniform site energy (i.e. with Sn = e Vra), but 
also for systems with local energy disorder. There, the varying coupling strength may 
effectively lift the disorder, since the localization depends on the ratio of detuning and 
coupling strength. 

It is also worth mentioning that in the regime where Jq ^ Is-, i-e. for a large 
decay rate into the sink, a quantum Zeno-type effect slows down, paradoxically, the 
absorption into the sink by projecting the system with high probability into a state of 
an unpopulated last site. Therefore, in this regime, variations of the inter-site coupling 
strength due to oscillations have no pronounced effect on energy transfer efficiency. On 
the other hand, in the regime of a weak sink rate and weak dissipation, the excitation 
undergoes many cycles before it is finally absorbed, thus amplifying the effect of the 
oscillations. 

In order to test, how a local detuning of the (at equilibrium position) uniform 
site energies affects the observed transport efficiency enhancement, we add an exciton- 
vibration coupling term, 

i7ex-vib = x{un+l{t) - Un{t)^ |^)(^|, (19) 

which also appears in a-helices [S3], to the Hamiltonian ([T]). The coupling strength 
X = 10 is chosen such that for a = 1/4 in the dimer, the energies are of the same 
order as the exciton coupling between the sites. Figure [7] (bottom,right) compares the 
transfer efficiencies with (dotted lines) and without this detuning (solid). Although the 
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Figure 9. Enhancement of energy transfer efficiency of the oscihating chain (breathing 
mode) as compared to the transfer efficiency of the resting chain coupled with Jq as 
both suffer strong dephasing at rate 7 = 1 with ai = N/\2 (left) and ai = 
(right), and 0i = 0. 



motion-induced efficiency gain suffers due to the detuning of the site energies and the 
thereby caused locahzation, the effect persists. 

The claim that the observed enhancement is indeed a quantum feature is tested 
by adding a dephasing environment, which in general leads to the loss of quantum 
coherence. This is modeled by an additional term in the master equation (|6]): 

TV 

LdepY^P = 7 J] (2^n O^n PO^n " {^n > P}) (20) 

n=l 

Since we expect quantum coherence to be a key requirement on which this effect relies, 
adding a dephasing environment should lead to a smaller enhancement. Indeed, the 
simulations in figure |8] show that the enhancement over the static case with Jmax 
disappears at a specific critical dephasing rate 7c. However, an enhancement of the 
transport efficiency over the static case with Jo still exists, even for large dephasing rates 
(see figurejo]). As expected, the critical dephasing rate becomes smaller as the system size 
increases due to the longer time over which the excitation is exposed to the coherence- 
destroying action of the environment. Please note, that in [191 EQ] an enhancement of 
the transport in the presence of dephasing is observed when adding dephasing to a static, 
disordered transport network. In contrast, the enhancement observed here comes from 
the concerted motion of the molecules, and it is an enhancement with respect to the 
static case with the same dephasing. This enhancement exists whether or not dephasing 
is present. 

4. Guided excitation transfer 

In the previous sections, we observed that motion and in particular oscillations can 
nontrivially interact with the quantum dynamics of the excitation such as to drive and 
enhance the energy transfer. We now extend the previous cases to an externally driven 
scenario. As a simple example, we investigate the sweeping of excited population guided 
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Figure 10. Left: Sink population for a Gaussian wave packet sweeping along the 
coupled chain of iV = 13 sites. The solid line marks the reference case without 
pulse, dashed/dotted lines indicate data for pulses of different speed, but same width 
<J — 1, and strength A — 1/6. Right: Transport efficiency gain due to the pulse when 
dephasing at rate 7 is taken into account. Curves are obtained for the optimal speed 
for each 7. 



by a Gaussian pulse that modifies the displacement and hence the interaction between 
the molecules as it travels along the chain. A toy model could be implemented in a 
simulation with trapped ions by shining a laser with a suitable spatial intensity profile 
onto the chain, for example. In a molecular context, a conceivable mean could be to 
attach the molecular chain to a nano-mechanical oscillator. 

We model the pulse such that the distance between molecules n and tt, + 1, 
dnit) = do — u(t)n,n+i, is changed from their equilibrium distance by 

^Ad„exp( J<"-^^r""' ). (21) 

where a is the pulse width, and v its velocity. The displacement due to the pulse remains 
Gaussian as the pulse uniformly propagates along the chain with its center at the position 
of site n, i.e. ndo = vt + do- The effect of this pulse is to concentrate the molecules 
in a certain region of the chain and thereby couple them more strongly, whereas the 
molecule density and coupling remains unaffected outside the pulse. Again, we would 
like to stress that there exist alternative handles to modulate the coupling strength, for 
example, a suitable change in the direction of the molecules' dipole moments. 

We start the pulse centered around site n = 1 at time t = 0, and first observe the 
moving Gaussian pulse according to (21) for different values of its speed v. All lengths 



are measured in units of do, such that we effectively set do = 1. In figure [lO] (left), we 
monitor the sink population as a function of time for different values of the pulse speed. 
The continuous line shows the reference case without a pulse, i.e. for the uniformly 
coupled chain with sites equally spaced at distance do- Dashed and dotted lines indicate 
solutions where a pulse is present. Clearly, there is an effective interplay between the 
moving Gaussian pulse and the dynamics of excitation transfer. 

When the speed of the Gaussian packet is zero, which means a stationary packet 
centered around the first site, it only compresses the distances in the beginning of the 
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chain (but leaves the remaining distances unchanged), and thereby causes stronger, but 
disordered couphngs. Although a stronger coupling increases the speed of the excitation 
transfer between two sites, for f = the efficiency of energy transfer, i.e. the sink 
population in the long-time limit, is less than that of a uniform chain without a pulse 
being present (solid line in figure 10 (left)). In fact, when the couplings at the beginning 
of the chain are stronger than between the remaining sites, they cause a localization of 
the excitation population at the beginning, due to the caused disorder in the couplings. 
Therefore, only a relatively small fraction of the population succeeds to reach the sink, 
whereas the localized population is finally dissipated into the environment. With this 
suppression, we once again face a distinct feature of coherent energy transport. In 
the classical, diffusive energy transfer, an increase of the hopping rates always causes 
an increase of the transfer rates, and consequently leads to a higher energy transfer 
efficiency. 

If the pulse moves, it will sweep along the excitation that is localized within the 
pulse. In the limit of very distant sites, a good strategy is to apply a strong pulse that 
couples only two sites and moves at the speed given by the Rabi frequency, such that 
after a full transfer to the second site, the pulse moves and couples site two and three. 
For a pulse that is too fast, only part of the excitation will be transferred to the next site, 
whereas if it is too slow, part of the excitation oscillates back onto the first site again. 
In the present case, we cannot neglect the coupling outside the pulse. Instead, we need 
to coordinate the two processes of localizing the excitation within the moving pulse, and 
the wave-like coherent transfer in the uniformly coupled chain. Figure 10 (left) shows 
that, for the given parameters and a chain length of 13 sites, the dominant contribution 
to the sink population due to the wave-like spreading of the excitation arrives between 
times 6-8 (in units of l/Jo)? with small additions at later times. The standing wave 
packet (dotted line) localizes the excitation at the beginning of the chain, but a fraction 
still succeeds to arrive in the sink at the same time as in the uniformly coupled chain. 
The increase of the sink population starts earlier than without the pulse being present, 
because the pulse also slightly increases the effective coupling strength for the rest of 
the chain, which leads to a faster propagation of the excitation. A wave packet moving 
with the speed of one site per time unit consequently arrives at the end of the chain 
after having passed 12 sites at time t = 12, where it leaves part of the excitation that 
it has swept along, resulting in an additional increase of the sink population. At the 
optimal speed of fopt = 2.53 ± 0.01, the wave packet arrives at about the same time 
as the excitation would in the uniformly coupled chain, hence effectively localizing the 
excitation in the front of the wave-like propagation. 

The optimum speed is proportional to the coupling strength, i.e. it increases linearly 
with Jo, because the Rabi frequency that is given by the coupling strength (up to units) 
is the speed at which an excitation moves between two sites during the Rabi oscillation. 
Therefore, the pulse has to match the speed given by Jq. 



The simulation is summarized in figure 11, which shows the increase of the 



sink population for a Gaussian pulse as compared to the uniformly coupled chain, 
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Figure 11. Energy transfer enhancement for a Gaussian pulse over the uniformly 
coupled chain with Jmax- Chain length is iV = 13, and A = 1/6. Contours of the 
surface (left) are plotted on the right. 



A = Psin\^{v,a) — -Psink(^max), i-G. at a fixed maximum intensity of the pulse. For 
small pulse widths and low velocities the gradient of coupling strengths changes quickly, 
and influences only few sites. Therefore, the efficiency suffers most due to the strong 
localization of the excitation within the pulse. On the other hand, for large values of a 
the wave packet is broad, that is, it covers a larger region of the chain, and the gradient 
of coupling strengths changes more smoothly, which leads to less details in the velocity 
dependence. In the limit of large pulse widths, the wave packet virtually covers the 
entire chain such that the reference case with Jmax is reached, and the efficiency gain 
approaches zero. We also found that there is a global optimum (wopt = 3.34 ± 0.01, 
Copt = 4.20 ± 0.05) where, at a given maximal intensity, the velocity and pulse width 
match best the coherent dynamics and achieve an increase of the transfer efficiency of 
more than 0.03 over the static chain with maximal coupling strength. Note that the 
local optimum from figure [lO] with a = 1 corresponds to the maximum of the front-left 



cut of the surface in figure 11 (left). 



Finally, we also account for a possible dephasing, as caused by fast thermal noise, 
for example, which is formally implemented by adding the term (20) to the Lindblad 



equation. The results are collected in the right panel of figure [l0| and show that even 
under substantial dephasing the gain due to the presence of the pulse persists, and 
may even attain values in the classically forbidden region of enhancements A > (not 
shown) . 



5. Concluding remarks 

In this work, we demonstrate general features of energy transfer efficiency in 
mechanically oscillating systems, and relate our model to the biological scenarios of 
coherent energy transport in proteins. Complementary to previous works, we focus 
on the driving of the excitation dynamics due to motion of the underlying molecular 
structure. The time dependence of the inter-site coupling arises from its distance 
dependence together with relative motion of the sites, for example. We generally 
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find a motion-induced quantum enhancement of tfie excitation transfer efficiency over 
tlie equilibrium configuration of tlie finear cfiain, and furtliermore over any static 
configuration tfiat is met during tlie mechanical motion. This effect distinguishes the 
quantum-coherent transport from the classical, diffusive transport. If biological systems 
manage to utilize the modulation of coherent couplings by their motion, they can profit 
from the resulting enhancement of the transfer efficiency. It would be interesting to 
extend our ideas beyond the linear setting to the enhancement of excitation transfer in 
light harvesting complexes in photosynthesis or light harvesting complexes in artificial 
devices. 
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